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Abstract 

The idea of glue clusters, i.e. short-range correlations in the quark- 
gluon plasma close to freeze-out, is used to estimate the width of 
balance functions in momentum space. A good agreement is found 
with the recent measurements of STAR collaboration for central Au — 
Au collisions. 



1 Introduction 

It was observed by Bass, Danielewicz and Pratt [T] that measurements of 
balance functions in nucleus- nucleus collisions provide a " clock" which allows 
to determine when the charges observed in the final state are created. The 
measurements performed by STAR O |3l H] and NA49 [5], |6] collaborations 
proved (or at least strongly suggested) that the charges appear at the late 
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stage of the production process, close to the freeze-out [7] , thus implying that 
the quark-gluon plasma phase is dominated by gluons. 

In addition, these experiments showed that the measured widths of the 
balance functions in central AA collisions are substantially smaller than those 
seen in pp and in peripheral AA collisions. The natural conclusion from this 
observation is that the presence of the quark-gluon plasma phase induces 
additional correlations in the system. The nature of these correlations is of 
course a matter of debate. They were extensively discussed [7] in the frame- 
work of the thermal blast wave model [8]. It was argued that, when the 
quark-gluon plasma is sufficiently cool, pairs of positive and negative charges 
are created, with thermal momentum distribution, in certain restricted do- 
mains in space regions (as implied by the HBT measurements [S]). The 
resulting momentum separation is measured through the balance function. 
Since the creation of charges happens close to the end of the evolution, the 
separation of charges induced by expansion of the system is not effective and 
thus a small width of the balance function is maintained. 

In [in], on the other hand, the additional correlations were interpreted 
as evidence for clustering in the quark-gluon plasma. It was shown that 
production of uncorrected isotropic clusters of gluons decaying into qq pairs, 
when supplemented by their coalescence into hadrons [11], explains -in a 
natural way- the correct width of the balance function in pseudorapidity. 

In the coalescence model [TT] it is assumed that hadrons are created by co- 
alescence of "constituent" quarks and antiquarks into qq, qqq and qqq colour 
singlets. If there are no correlations in the system before coalescence, then 
the only correlations which can appear between hadrons can be those induced 
by resonances (which also result from coalescence). This picture describes 
reasonably well the hadron-hadron data [12J. The additional short-range 
correlations, needed to account for the experimental experimental measure- 
ments of STAR [21 [3l 15 and NA49 [3 [6] collaborations, may only appear if 
the constituent quarks and antiquarks are correlated between themselves be- 
fore the coalescence process takes place. The postulated glue clusters supply 
these additional correlations. 

Since pseudorapidity is determined solely by the particle production an- 
gle, distributions in pseudorapidity are only marginally sensitive to particle 
momenta. It is therefore interesting to verify if the hypothesis of glue clusters 
can explain also the recently measured [3] balance functions in momentum 
space. This is the subject of the present investigation. Our main conclu- 
sion is that, indeed, the presence of glue clusters at the last stages of the 
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evolution of the quark-gluon plasma can account for the observed width of 
balance funtions in qside, Qiong and Qout- This result requires that the average 
momentum < q > carried by a quark and antiquark in the decay of a glue 
cluster is located around ~ 120 MeV. 

In the next section we derive the momentum dependence of the balance 
function. Comparison with data is presented in Section 3. Our conclusions 
and comments are listed in the last section. The Appendix summarizes our 
treatment of acceptance corrections. 



2 Balance function from glue clusters 

Consider a two-body decay of a glue cluster into a qq pair. The distribution 
of the decay products is 

djTi 

Vc{Pl-V2)KVl+V2-P) (1) 



dpidp: 



where P is the momentum of the cluster, pi and p2 are momenta of the decay 
products and the function Vc is responsible for the details of the decay {dp 
stands for three-dimensional relativistic phase-space, dp = d^p/E = d'^p±dy). 

Consequently, the distribution of quarks and antiquarks arising from two 
clusters is 



dN{pu,Pu,Pd,Pd) 



dPudPDPc{Pu)Pc{PD) 



dpudpadpddpa 

S{pu +Pu- Pu)vc{Pu - Pu)S{pd + Pd - PD)Vc{Pd - Pd) (2) 

where pc is the distribution of clusters in momentum space. The subscripts 
f/, D refer to the flavour of quarks to which the cluster decays, and m, c?, u, d 
to flavour of the quarks and antiquarks themselves. 
For the balance function we obtaiij^ 

B{j>+,p^) j dpudpudpddpd^ipu + Pd - P+)KPd + Pu - P-) 

.dN{pu,Pu,Pd,Pd) .ON 
^{Pu-Pd)^{Pd-Pu) , , , , 3 

dpudpudpddpd 

where describes the coalescence process. 



"'^ Since we are only interested in the width of the balance function, henceforth we shall 
ignore all normalizing factors. 
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Introducing 

p+ = Pu + Pd ; p- = Pd + Pu] 6+=pu- Pd ; S-=Pd-Pu (4) 

we have 

B{p+,p-) ~ I d5+d5^p,{pu + Pu)pciPd + Pd)^{S+)^{S.) 

Vc{Pu - Pu)vc{Pd - Pd)S{Pu + Pd -P+)S{Pd+Pu-P-) (5) 

and thus 

5(5) ~ J dS+d5-^5+ + 5-)/2]<l>[5+ - 5-)/ 2] 
J dpp^iip + r )/2]p,[(p - 6-)/2]v,[{6 + 6^)/2]v,[{6 - 5+)/2] (6) 

where 

P = P++P^; 6 = p+-p^; 6^ = 5+ ±6^. (7) 

Eq. ([6]) represents our final resuh. Its consequences for the widths of the 
balance functions in various components of momentum are discussed in the 
next section. 



3 The widths of balance functions 

To simplify the discussion and reduce the number of parameters, we assume 
in this section that the functions Vc{x) and Pc{x) are Gaussians: 

v,{x) ~ e-^"/^^ p,{x) = e-i/i-^i/^l (8) 
Then the dependence on and 6~ in ([6]) factorizes out and we have 

B{5) ~ e-'" n{Sy, n{S) = J dpe-^V'-le-Py^'" (9) 

If there are no acceptance restrictions, the integral extends over full phase- 
space, Vt{5) is a constant, and the balance function is isotropic: its width does 
not depend on the chosen directioiJl and is entirely determined by the pa- 
rameter v'^. In the STAR experiment [3], however, the acceptance corrections 
are important. They were estimated and are summarized in the Appendix. 



^This was already observed in [7]. 
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Eq.(|n]) contains three parameters: f, r and ry. This number can be still 
reduced, using the available information on the distribution of transverse 
and longitudinal momenta. Indeed, the single particle transverse momentum 
distribution can be expressed as 

D{p+^) = e-2p+x/('-'+-')rfV± ^< pI± >= (r' + v^)/2 (10) 

with an analogous formula for pj| . Consequently we have 

r'^ = 2 <pI> -v^; rj = 2 < > -v^ (11) 

and we are left with one parameter, f ^ which should be adjusted to data. 

One sees from these formulae that ry is very large and therefore its exact 
value is not important for estimate of Siong- Therefore, to determine v, we 
used the value < 6iong >= 190 MeV, measured in the central Au — Au 
collisions [3]. Using Eq. ( 123|) from the Appendix, one obtains v = 276 MeV. 
With this value and < p± >= 400 MeV [E!, Eqs. (EI]) and ([18]) give 

< Ss^de >= 284 MeV; < 5out) >= 126 MeV (12) 

This should be compared with 

< 6side >= 280 ± 10 MeV; < Sout) >= 0.110 ± 10 MeV, (13) 

measured in [3]. Given the crudeness of our (Gaussian) approximation, the 
agreement is more than satisfactory. The width uncorrected for acceptance 
is < 5 >= 220 MeV. 

From the obtained value v = 276 MeV and using ([T]) one can evaluate the 
average value < g > of the momentum carried by the quark and antiquark 
in the decay of the cluster. One obtains < q >= v^jA: ^ 122 MeV and 
v/< > = v/2 ^ 138 MeV. 

4 Conclusions and comments 

In conclusion, we verified that the hypothesis of glue clusters, i.e. positive 
short-range correlations in the quark-gluon plasma [10] can account for the 
small widths of the balance functions in momentum space, observed recently 
by the STAR collaboration j3] for central Au — Au collisions. In particular, 
it was shown that when the parameters of the model are determined from 
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the observed < Siong >, one obtains reasonable values of < Sgide > and of 
< ^out >■ This result confirms the existence of correlations in the plasma. 
It also shows that they can be effectively studied by the method of balance 
functions [Ij. 

Several comments are in order. 

(i) It should be emphasized that the observed small width of balance func- 
tions in momentum space simply implies the existence of additional short 
range correlations between particles produced in central nucleus-nucleus col- 
lisions (as compared to pp collisions). The nature and origin of these corre- 
lations remains a subject of debate but their very existence is beyond doubt. 
For example, the explanation of the balance functions in the blast wave model 
[H], presented in [7], also exploits the correlations in the plasma (in the form 
of " domains" from which the balancing charges are emitted) . 

(ii) The existence of short-range correlations in the quark-gluon plasma 
close to the phase transition should not be surprizing. Indeed, the lattice 
calculations indicate strong deviations from the Stefan-Boltzman limit even 
at temperatures largely exceeding Tc [H] . Thus the presence of quasiparticles 
is likely. The nature of these quasiparticles is an open question. Our results 
suggest that they may take the form of gluonic clusters. 

(iii) Although our semi-analytic estimates are admittedly rather crude, it 
is interesting to observe that the average momentum characterizing the decay 
of a glue cluster (120 MeV) seems close to the break-up temperature in the 
blast wave model, as discussed in [7]. This may suggest that the approach 
of [7] and our interpretation may not impossible to reconcile. Precise data 
on balance functions for KK and Ktt pairs may throw some light on this 
problem. 

(iv) In our argument we have used the coalescence model [TT]. Although 
the model is rather successful in explaning, e.g., the particle content [15] 
and V2 scaling [16], its basic idea is often questioned since it is difficult to 
reconcile with the picture of pion as the Goldstone boson of chiral symmetry 
breaking [17]. A possible way out is to admit that the constituent quark 
mass depends on the temperature of the system (or rather on its distance 
from the temperature the chiral transition). If true, then also the mass of the 
glue clusters will depend on the temperature. Although our results do not 
depend on the constituent quark mass, it may be interesting to investigate 
this problem in more detail. 

(v) It was noticed in [7] that corrections due to production of hadronic 
resonances may change the results by 10-20 %. This is not dramatic for 



6 



our semi-quantitative approach but must be eventually improved if serious 
comparison with data is attempted. 

It is well known that resonance production is essential in description of 
the short range correlations in pp collisions [12] and that it allows to explain 
the single particle spectra in nucleus- nucleus collisions (see e.g. [18l|)- Fur- 
thermore, it has been shown [TH] that thermal, uncorrelated production of 
resonances cannot reproduce the small width of balance functions in central 
nucleus-nucleus collisions. It should be noticed, however, that if hadronic 
resonances are formed by coalescence of q and q which are decay products of 
glue clusters, they are correlated. This -in turn- should produce additional 
correlations between their decay products and thus reduce the width of the 
balance functions. It would be very interesting to investigate such a possibil- 
ity in detail. This, however, requires much more sophisticated analysis than 
the one presented here. 

(vi) The present work represents only a first order approximation to the 
problem. We only considered the simplest case when the hadron pairs are 
created by coalescence of quarks and antiquarks which are decay products 
of just two clusters. The contributions from more than two clusters should 
certainly be included in a more precise analysis, although it seems likely 
that they are suppressed because (i) the probability to have more clusters 
very close in space (as required by the coalescence mechanism) is small and 
(ii) matching the colour of quarks from two different clusters reduces this 
contribution even further. At present we can only mention that (a) for three 
clusters the width of balance function depends on details of the coalescence 
process and therefore its estimate involves more parameters than hitherto 
considered, and (b) if clusters are uncorrelated, contributions from 4 clusters 
cancel in the balance function, see e.g. [TU]. 

Let us also add that the Gaussian approximation is certainly rather crude 
and our treatment of the acceptance corrections is at best approximate. This 
can of course be improved, although we feel that it would be not justified at 
the present level of understanding. 

In summary, the present investigation shows that the hypothesis of glue 
clusters in quark-gluon plasma can account for the recently measured widths 
of the balance functions in central Au — Au collisions. This may have im- 
portant consequences for the phenomenology of the quark-gluon plasma. It 
remains an open question if this result is confirmed when data in larger ac- 
ceptance regions are available. 
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5 Appendix: Acceptance corrections 



We first consider transverse directions, i.e. two-dimensional vectors p+ and 
p_. The distribution is given by (Q with the acceptance hmits 

< pI. + Ply < A^; < p\ + p\ < A' (14) 

where D= 200 MeV and A= 600 MeV. 

We are interested in the distribution of Sout and Sside, defined with respect 
to total momentum of the pair. Since the system is invariant with respect to 
common rotation of all vectors, one can select the y axis as "out" and the x 
axis as " side" . Thus we are looking for the configuration in which 

P^+ +p^- = Sside = 4 = Px+ - Px~ = '^Px+; P±y > 0. (15) 

Consider first 6 out = Sy. This is the most complicated case because in 
this direction there is a significant effect of the flow which strongly distorts 
the spectrum and must be taken into account. Flow implies that the cluster 
moves (on average) with the velocity V of the fluid. Therefore the labo- 
ratory system in which it is observed moves with the velocity —V and the 
distribution of the decay products is 

e-(^L.)V2.^e-(^^-)V2^'^ C = - P'-; Pout = Pi + P'^ (16) 

where 

pI = i\p±-VE±]- E± = ^pl + x^ + ml (17) 

Consequently, the observed distribution of the balance function is (we abbre- 
viate p^x = 

P{5out) = r dxe-'^'/^' e-^-/2''^e-^"-/2'^^ (18) 

JO Pmin 

where Pmax = 2A{x) - 6out, Pmin = 2D{x) - 6out for 6out < 2-D(x), Pmin = 

for 6out > 2D{x) and 

A(x) = VA2 - x2; D{x) = VD^ - x2 (19) 

Numerical evaluation of ( fTSl) shows that, with V = 0.6, the effect of flow 
reduces significantly the width. 
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For Sside = the effect of flow is largely cancelled and we neglect it. 
Consequently, the distribution is 



Jd{x) 



^ g -{P+y+P-yr/2r^f,-{P+y-P-yr/'^^'' 



(20) 



where x = 6 side f^- One gaussian integration can be expressed in terms of 
error functions and one obtains 



Jd(x] 



erf{s\B+\) + erf{s\B_\) 



A A_ 
-r^er/(s|^+|) - er/(s|^_|) 



'■+1 



(21) 



with 



Aj^^±D{x) +w; S± = ±A(x)+«;; w ^ y+ 



T — V 

J.2 _|_ y2 



r + V 



,(22) 



Note that we always have > and B_ < 0. 

Finally, let us consider 6iong- In this case there is no lower limit on particle 
momenta. Consequently, the distribution becomes 



P{Slong) ~ e 



Jo 



(23) 
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